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The — > L>*^z/ form factor at zero recoil and the determination of |V^,| 
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We summarize our lattice QCD study of the form factor at zero recoil in the decay B — > D*lu. After careful 
consideration of all sources of systematic uncertainty, we find, /ia^I) = 0.913^*l3Q, where the first uncertainty 
is from statistics and fitting while the second combined uncertainty is from all other systematic effects. 



1. INTRODUCTION 

A precise value for CKM matrix element \V c b\ 
is an important ingredient in the study of CP vi- 
olation and the determination of Wolfenstein pa- 
rameters (p, fj) via the unitarity triangle. 

Experimental studies of the B — > D*iv decay 
rate determine the combination \V c b\ /^(l)!]]] [§■ 
The hadronic form factor at zero recoil 

must be computed by theoretical means in order 
to extract \V c b\. 

Heavy quark symmetry imposes powerful con- 
straints on /i J 4 1 (l): it requires /^(l) — > 1 in the 
infinite mass limit Q and it determines the struc- 
ture of 1/wq corrections^. These power correc- 
tions correspond to long distance matrix elements 
in Heavy Quark Effective Theory (HQET). Previ- 
ous determinations of /^(l) have employed sum 
rules H or have appealed to quark models |J. 

We present a determination of h Al (1) from lat- 
tice QCD using a new double ratio method de- 
signed so that the bulk of correlated statistical 
and systematic errors cancel M. In fact, with this 
method nearly all errors - including quenching - 
scale with 1 — h Al rather than with h Al ■ 

Our result is computed in the quenched approx- 
imation. The method lends itself to a completely 
model independent determination of h Al {\) once 
unquenched gauge configurations become avail- 
able. The full details of our analysis are found in 
Reference [pi. 



2. THE METHOD 

In the heavy quark expansion! 



h Al {\) = j] A 



1 



,(3) 
-v 



(2m c ) 2 (2m c )3 {2m b f 



• p 



2l A 



( 1 



1 



(2mb) 3 2m c 2mb 2m c 2mb \2m c 2mb 



»(3) 



+ 



/ 1 



1 



2m c 2m,b \2m c 2mi 



+ 



(1) 



The short distance quantity n A relates HQET to 
QCD and is known to two loop order |1. The 
£'s are the matrix elements in HQET which we 
determine from lattice QCD. 

We construct double ratios Q from the lattice 
matrix elements (D | V 4 | B), (D* | A x \ B) 
and (D* | V4 | B* ). The matching between lat- 
tice and HQET currents is known to one loop 
order By varying the "charm" and "bottom" 
quark masses in the lattice simulations, we are 

able to extract all of the required Cs in Eq. Ill, 

(3) . . . 

except £_ . The truncation uncertainty in our 

h Al (l) determination is the size of this unknown 

term plus higher order terms in the heavy quark 

expansion. 

3. RESULT FOR h Al (l) 

We find 
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h Al (1)=0.913 
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Tabic 1 

Budget of statistical and systematic uncertainties 
for /^(l) and 1 — /^(l). The row labeled "to- 
tal systematic" does not include uncertainty from 
fitting, which is lumped with the statistical error. 
The statistical error is that after chiral extrapo- 
lation. 



uncertainty 




(1) 


1 - h A 

% 


(1) 


stats and fits 


+0.0238 


-0.0173 


+27 - 


20 


adj. m c & mb 


+0.0066 


-0.0068 


+ 8- 


8 




±0.0082 


± 9 




a s {K/2m Q ) 2 


±0.0114 


±13 


(A) 3 /(2m Q ) 3 


±0.0017 


± 2 




a dependence 


+0.0032 


-0.0141 


+ 4- 


16 


chiral extrap. 


+0.0000 


-0.0163 


+ 0- 


19 


quenching 


+0.0061 


-0.0143 


+ 7- 


16 


total syst. 


+0.0171 


-0.0302 


+20- 


35 


stat © syst 


+0.0293 


-0.0349 


+34- 


40 



The first uncertainty is the statistical error added 
in quadrature to the uncertainty from fitting pro- 
cedures. The second uncertainty is the combined 
error, in quadrature, from the other sources listed 
in Tab. [jj, 

3.1. Systematic Uncertainties 

• Statistics and fitting procedures. The sta- 
tistical error in our result is after the "chiral" 
extrapolation of the spectator quark. The fit- 
ting procedure error includes the effect of excited 
state contamination and also bounds variations 
in /i^i 1 (l) from alternate plausible time ranges 
in fitting three-point ratio plateaus. Minimum 
Chi-square fits were obtained for all fits. Fits 
include the data correlation matrix and produce 
bootstrap error determinations. Small poorly de- 
termined eigenvalues of correlation matrices were 
rejected in a SVD decomposition. 

• Adjustment of m c and mb- The charm 
and bottom masses are determined by adjusting 
the bare quark mass until the kinetic mass of a 
lattice meson matches the physical D s and B s 
masses respectively These kinetic masses tended 
to be quite noisy. Our uncertainties are taken to 
encompass our m c and mb determinations from 
quarkonia. The agreement between heavy-light 



and -onia quark mass determinations is better 
for charm than bottom. Hence, our charm quark 
mass has the smaller uncertainty. 

• Matching, a^. Although the matching be- 
tween HQET and QCD is known to two loop 
orderj9|, the matching between lattice and QCD 
is only known to one loop order flo[| . Hence, we 
match among HQET, lattice and QCD schemes 
consistently to one loop order, choosing scales for 
the QCD coupling according to the BLM pre- 
scription The uncertainty in Tab. [I] reflects 
our estimate of non-BLM terms at and beyond 
two loop order in the perturbative expansions. 

• Undetermined (A/(2mq)) 3 terms. We did 

not determine £_ in Eq. [j]. Taking nominal val- 

iks for m c and mj and £_ « A rs 0.5 GeV 
we estimate the size of the unknown term to be 
~ 0.0017. This estimate is consistent with the 
size of all other 1 /rn,Q terms that are included in 
our result. 

• Action and currents a s (A/(2mQ)) 2 . The 
action and currents are tuned to tadpole- 
improved tree level. A careful analysis(|] of the 
form factors shows that remaining corrections are 
of order a s (A/(2m,Q)) 2 . Using nominal values, 
we estimate a s (A/(2m c )) 2 « 0.008. We also es- 
timate this uncertainty by repeating our analy- 
sis using tree level quark masses instead of the 
"quasi one-loop" masses used in our standard 
analysis. Wc estimate the uncertainty in /i^ 1 (l) 
to be ±0.0114 by this method. 




fl lpls (GeV ) 
Figure 1. Lattice spacing dependence for /^(l). 
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• Lattice spacing dependence. Wc computed 
/lAi(l) for a strange-mass spectator quark at 
three lattice spacings corresponding to (3 = 6.1, 
5.9 and 5.7. T he la ttice spacing dependence is 
shown in Figure 3.1. We take the weighted aver- 
age (horizontal line and error envelope) of values 
at our two finest lattice spacings as our best deter- 
mination of \ia x (1) ■ The value from our coarsest 
lattice is not used in this average since it from suf- 
fers larger heavy quark discretization errors. We 
do include the coarsest lattice in our bound on 
discretization errors. This bound is shown by the 
linear fit to all three points. 

• Chiral extrapolation. We have studied the 
spectator mass dependence for at our two 
finer lattice spacings. Figure |3~l| shows the linear 
extrapolation to the physical down quark mass 
for (3 = 5.9. We observe a similar slope, in di- 
mensionful units, for /3 = 6.1. Our best value 
of 1ia 1 (1) with a strange-mass spectator quark is 
shifted downwards by the amount shown in the 
figure to yield our final value of /i J 4 1 (l) with a 
bottom spectator quark. Statistical errors are in- 
creased as indicated in the figure. 
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Figure 2. Spectator quark mass dependence at 
(3 = 5.9. 



Randall and WiseJT| have computed pion loop 
effects upon h,A 1 (l) in Chiral Perturbation The- 
ory. The curve with the cusp shown in Fig. [3.l| 
is this prediction for the spectator mass depen- 



dence^. The departure from linear behaviour 
near the down quark mass adds an additional un- 
certainty to the chiral extrapolation. 
• Quenched QCD. Our result is obtained in 
quenched QCD. The quenched coupling runs in- 
correctly: the short-distance quantity t\a changes 
by +0.0050 when quenched. Long distance form 
factors, such as hA 1 /r/A, are also affected by 
quenching. With our method, however, quench- 
ing errors only affect the deviations from unity. 
Guided by studies of decay constants fs and /d, 
we expect this error to be less than 10%. Our 
quenching uncertainty reflects estimates of both 
long and short distance effects. 

4. |V cb | 

Using our result, wc find: 

45.9±2.4±^ CLEO[l] 
38.7 ± 1.8±1;| LEP§ , (2) 
39.3 ± 2.5t^ Belle [|] 



10 3 |V f 



where the second error results from adding all our 
uncertainties in quadrature. This result includes 
a +0.007 QED correction to /^(l). 
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